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1. Introduction 

In the context of current efforts to build quantum networks [1], a promising way to 
transfer quantum states reliably in the network is through the coupling of single photons 
and atoms in the setting of cavity quantum electrodynamics (QED) [2]. For practical 
applications, the coupling between a single atom located in an optical cavity and a 
single intercavity photon should be strong. The strong coupling condition requires that 
go/2 ^ r, K, where go is the one-photon Rabi frequency, T is the atomic decay rate to 
modes other than the cavity mode, and k, is the decay rate of the cavity mode itself. 
This may be achieved using high-finesse optical cavities, with an extreme reduction in 
the cavity volume, and using atomic transitions with a large dipole moment. 

Using the cavity QED techniques, schemes for a number of operations of direct 
relevance to quantum information processing have been proposed, one of them being 
generation of single photons 'on demand'. Through strong coupling of a mode of the 
cavity field to an atomic transition, which is resonantly driven by the input single- 
photon pulse, an external control field of Rabi frequency ^Icit) transfers one photon in 
the cavity mode to a long-lived atomic memory, which can then be released at will to 
free space through the cavity output mirror, leading to an output single-photon pulse as 
a collimated beam. The temporal structure (both amplitude and phase) of the resultant 
'flying photon' can be tailored by way of the control field ficit) [31 H], with the spatial 
structure of the wave-packet being set by the cavity mode. 

The basic scheme (see figure 1) involves a three- level atom in a A-configuration 
with an excited state \a) and two lower states \b) and |c). An optical cavity mode is 
strongly and coherently coupled to the atom on the \b) ^ \a) transition with rate go, 
and a strong classical field fic(t) drives the |c) \a) transition. Denoting by \x,n) 
a state in which the atom is in state \x) and there are n photons in the cavity mode, 
reversible transfer of a state between light and a single trapped atom can be achieved 
through the mappings \b,l) -H- |c, 0) for the coherent absorption and emission of single 
photons by a procedure involving the 'dark state'. The atom-cavity system coupled to 
a classical control field Qc{t) has an instantaneous eigenstate \D), given by 

\D) =cos9{t)\b,l) -sme{t)\c,0), (1) 

where 

ta„.W = ^. (2) 

([1]) is called the 'dark state' [5] as it does not have any component involving the excited 
level |a), which can radiate. If initially the system is prepared in the state |c, 0) with 
Qc(t = 0) = 0, then a sufficiently slow increase of Oc(t — > oo) to Qo ^ go leads the 
state to move adiabatically to \b,l). On the other hand, if initially the system is in state 
\b, 1) with flcit = 0) = f2o, then decreasing flcit) adiabatically to zero leads the system 
to the state |c, 0). The advantages of the dark state protocol are: (a) it minimizes the 
dissipative effects generic to two-level systems, and (b) it is reversible, i.e., a photon that 
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Figure 1. Three- level A scheme 



is emitted from a system A is efficiently transferred to another system B by applying 
the time-reversed (and suitably delayed) control field flc(t) to system B. 

Cirac et al. [3] proposed a quantum network with nodes consisting of cavities 
each containing a three-level atom. The entanglement between atom-cavity states 
at each node is generated using the protocol described above. Similar protocols 
involving adiabatic transfers using three-level atoms, either trapped in or passing 
through electromagnetic cavities, have also been proposed for other purposes. Parkins 
et al. [Ej were the first to propose this kind of protocol to generate Fock states and other 
nonclassical states of the cavity mode. Pellizzari et al. [7] proposed implementation of 
a two-bit quantum gate by putting two three-level atoms inside the cavity. The first 
experimental attempt of such a reversible mapping of a field to and from an atomic state 
has been made [8] by using a single trapped cesium atom. In this example, |c) and \b) 
represent internal states of the atom with long-lived coherence, namely, the hyperfine 
states in the 6S1/2, F = 3 and F = 4 manifolds of atomic Cs, and \a) corresponds to 
6P3/2, F = 3. Intrinsically reversible and controlled single-photon sources have been 
demonstrated P [TOl [TTl [12] using the dark state based on stimulated Raman adiabatic 
passage between two ground states of a single atom strongly coupled to a single mode 
of a high-finesse optical cavity. 

The above kind of adiabatic transfer has also been used with ensembles of 
three-level atoms in a pencil geometry, which lead to the well known phenomena of 
electromagnetically induced transparency (KIT) and slow light [T3 | [T^ [T5l [T6| [T7 t [T8l [T9] . 
With atomic ensembles, the light pulses can be stored as collective atomic polaritons 
and recovered with high fidelity again by manipulating a control laser adiabatically [20] . 
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Duan et al. [21] have proposed long-distance quantum communication and entanglement 
using nodes of atomic ensembles. Very recently, single-atom EIT condition through a 
coherent dark state has been achieved in a high-finesse optical cavity [221 ES]. The 
single atom effectively acts as a quantum optical transistor, coherently controlling the 
transmission of light through the cavity. 

Though the various procedures outlined above have different physical features, 
underlying all of them is the assumption of adiabatic transfer through the dark state 
of a single three-level atom in the A-configuration. The adiabatic condition is taken to 
be fulfilled if the evolution time is significantly longer than the inverse of the frequency 
gap between the dark state and the other eigenstates of the Hamiltonian. But it is of 
utmost interest to understand the precise conditions, the fidelity and the experimental 
limitations of this procedure which allows transfer of quantum state information from 
the photon to an atom and vice versa. 

This question has received a lot of attention in several physical contexts. In 
general, the superadiabatic transfer procedure given by Berry [21] allows one to treat 
the non-adiabatic corrections in the quantum evolution for slow rates of change of the 
control laser in the absence of dissipation. For cavity-QED with a trapped atom, Duan, 
Kuzmich and Kimble [25] have presented a detailed analysis of schemes which require 
photon transfer between cavity mode and the mode of the external channel. Here the 
non-adiabatic analysis requires inclusion of a large number of states, and the analysis 
has been carried out numerically. A drawback of this analysis is that the dissipative 
terms are put in the Schrodinger equation for the amplitudes. This may suppress some 
important non-adiabatic elements of the quantum transfer process, and the issue needs 
to be examined, if possible, using a better formalism. The same problem has also been 
analyzed by Yao et al. [26], but in their analysis, the coupling to the external channel 
has been treated in the Wigner-Weisskopf approximation, which may not hold for all 
the procedures of interest. 

We wish to address here the exact problem of the transfer dynamics in the context 
of a single atom in an electromagnetic cavity, for arbitrary rates of variation of the 
control laser and with the inclusion of dissipation in the form of spontaneous emission 
from the upper level. For a three-level atom in the A configuration interacting with a 
single cavity mode, one can carry out a thorough non-adiabatic analysis as the Hilbert 
space is small. As has been seen in [^, the effect of spontaneous photons from the 
excited atomic state on quantum computation is much more destructive than cavity 
decay. This is because after a cavity photon emission, the system is still in a dark state, 
which is not necessarily true for spontaneous emission. Thus the atomic excited states 
can be populated, which gives rise to further spontaneous photons. Furthermore, we 
work in the strong-coupling regime where Qq 2> k. In this regime we neglect the decay 
of the cavity photon. Thus our analysis is not directly applicable to cases where photon 
transfer to external channel is of comparable time scale. 

The paper is organized as follows. In section 2, we first study the problem of an 
atomic system in a A configuration interacting with photons in a cavity mode and a 
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time-dependent control field, ignoring any dissipative effects. Our simple considerations 
show that the non-adiabatic effects are large when Qc{t) is small. Since during storage 
as well as retrieval, flc{t) is made zero, the non-adiabatic effects are unavoidable, and 
the problem cannot be treated perturbatively. Accordingly, we study the problem 
numerically. 

Next, in section 3, we analyze the effect of dissipation by allowing for the possibility 
of spontaneous emission, which is inevitable in any realistic system. Since spontaneous 
decay sends the system to the lower levels, one might expect that dissipation would 
mitigate the non-adiabatic effects. We study the impact of dissipation using a method 
in which the quantum evolution of the wavefunction is interrupted by spontaneous 
transitions of the state to the lower levels \b) and |c). Spontaneous decays are governed 
by a stochastic Poissonian process with a rate which we take to be the width of 
the level \a). Though dissipation in quantum systems is naturally incorporated in 
the density matrix formalism, we believe that the wavefunction treatment is a fair 
approximation, and it has two distinct advantages. First, the physical picture of the 
system evolution and the role of dissipation is rather transparent. Second, there is a 
numerical simplification of having to solve only three coupled time-dependent differential 
equations. 

In section 4, we present our results on the dynamics of the fidelity of the storage and 
the retrieval process. Different rates of variation of the control field at different signal 
strengths are explored, each in the absence and presence of dissipation. We compare 
our results in the absence of dissipation to the superadiabatic theory due to Berry 
[24] . Finally, in section 5, we present our conclusions. In Appendix A, we compare the 
wavefunction approach with the density matrix approach and point out the nature of 
the approximation made in our analysis. 



2. Formulation with a single isolated atom 



We now consider the A atomic system as shown in figure 1, which interacts with the 
signal photon in the cavity mode and the control laser field. The Hamiltonian of the 
isolated atom-signal field system is Hq + Hj, with 

1- 



Ho = a^a^ + - + XI Ex\x){x 



(3) 



where a^ denotes the annihilation operator for the particular cavity mode of frequency 
u, which couples to the transition between levels \a) and \b), and x = a,b,c. In the 
rotating-wave approximation, the time-dependent interaction Hamiltonian of the fields 
with the atomic system (in one dimension) is 



Hi = %o \a){b\ais + aLl\b){a 



^cWe |c)(a 



h 



nc{t)e-"''\a){c 



(4) 
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where v is the frequency of the control field, and (70 ? the Rabi frequency for the cavity 
photon is given by 



with ^ab being the dipole moment between levels \a) and Vq the cavity volume, 
and eo the permittivity of vacuum. It is assumed that only transitions \a) — )■ \h) and 
I a) — )■ |c) are dipole-allowed. The control field with a Rabi frequency VLq is treated 
classically. 

Considering an n-photon quantum field, the wavefunction of the system of one 
atom+field can be expressed in general as 

\^>{t)) = J2 [An{t)e~'''-'\a,n) + Bn{t)e-'''>'-'\b,n) 

n 

+ C„(t)e-*^-*|c,n)]. (6) 

Here 



2, 

For this one-atom closed set of equations of motion for the coefficients An, Bn 

and Cn are obtained. For further analysis, it is convenient to work with the vector X 
with components 

ttn = v4„, bn = G ^ Bn+1, Cn = G Cn, (8) 

where As = u — Uah and Ac = u — Uac denote, respectively, the detunings of the two 
fields. Then X obeys the equation 

= H(t) X(t), (9) 



dt 
where 

/ Qn -^cit) \ 

H(t)= Qn As 

\-nUt) Ac I 



(10) 



with Qn = go\/n + 1 denoting the Rabi frequency of the n-photon signal field. 

We write the solution of the above equation in terms of the instantaneous 
eigenvalues and eigenvectors: 

H(t) = n\k{t) \uk{t)), k = l,2, 3. (11) 

At two-photon resonance, i.e., with As = Ac = A, the instantaneous eigenvalues are 

Ai = A, A2 = y + f]R(t), Xs = j-nR{t), (12) 

where 



n,K{t) = J9l + \nc{t)\'. (13) 
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The instantaneous eigenvectors of the atom-field system are 
\uiit)) = cose{t)e'^ 1) + sin^(t) |c,n), 

\u2{t)) = COS ^^^^ |a, n) + sin^^^ sin6'(t) |6, n + 



-cos^(t)e~** \c,n) 

m 



sm 



\a,n) + cos 



mi 



sm9{t) \b,n + 1) 



-cos^(t)e-*^ |c,n)". 
Here (j) is some arbitrary constant phase of the control field, and 
tane{t) 



9n 



A/2 ' 



(14) 

(15) 
(16) 



The eigenstate \ui{t)) is the 'dark state'. Now we can expand the solution in terms of 
these eigenvectors as 



k 

The time-dependent coefficients D^it) obey the following equation: 
dD^it) 



(17) 



dt 



+ Dm{t){Um{t)\Um{t)) 



= - ^D,(t)(n^(t)|u,(t))e-^/o(^^(*')-^-(*'))'^*'. 

By making a change of variable to 

VUt) = DUt)e'IoM''^''', iP.,{t') = (n^(t)|n„.(t)), 
the evolution of Vmif) is obtained as 



dVm{t) 

dt 



- E V^fc(t)(^n^Wl^fc(t))e-^/o^-(*')^*', 



:i8) 



(19) 



(20) 



where A'^^ = Xkm + Pkm, hm = Afc-A,„ and Pkm = Pk-^m- Note that the assumption of 
adiabaticity implies that the coefficients DmS or VmS are independent of time. The time- 
variation of these coefficients is governed by the terms {um{t)\uk{t)) , which essentially 
give rise to non-adiabatic effects. From the eigenvectors in Eqs. (fT4l) . we derive the 
following: 



Ml M2) = 


- {U2\ui)* 


= 6 it) sin 


Ml M3) = 


- {uslui}* 


= e{t) cos 


M2 M3) = 


- {U3\U2)* 


m 

2 ' 



-icj) 



(21) 
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where 



m = — 2 (22) 



^ -it- P^' 

All /3/;S are zero, and therefore A'^^ = \km- 

Then, finally we can express the coefficients C„(t) of our initial 

basis of the bare states in ^ in terms of these solutions as 

A^ft) =(cos^\/2(t)e-^/o^'.(*')'^t' 

-sin^\/3(t)e-'/o^3(*')'^*'), (24) 
5„+i(t) = yi(t)e-'/o cos^(t)e''^ + sin^(t) 

X (sin^\/2(t)e-^/.^^(*')^*' 

+ cos^r3(t)e-/o^3(0^*')^ (25) 

C„(t) = \/i(t)e-' /o sin e{t) - cose{t) 

xe-^(sin^\/2(t)e-^/o^2(0<^*' 

+ cos^V,it)e-^I>'^^''^'''). (26) 

From the above equations, it is quite evident that the non-adiabatic perturbation, 
which is proportional to 6, becomes large in the storage and retrieval process as flcit) 
becomes small. Thus we solve Eqs. (12 Op for VmS numerically. For simplicity, we take the 
phase of the control field, = 0, and the optical detuning, A = 0. The latter would 
imply that = f and ^Z' = 0, leading to considerable numerical simplification. 

It is interesting to place the significance of our analysis in the context of a known 
general result on adiabaticity for three- level systems obtained by Oreg et al. [27] . These 
authors analyzed the density- matrix equations as SU(3) rotations of an 8-component 
vector S constructed out of 8 independent components of the density matrix. It was 
shown that the stationary sector of the solution consists of two vectors Fi and F2, 
which are obtained from the Hamiltonian parameters. This implies that 5 = 0, when S 
is any linear combination of Fi and F2. Such an S follows the subspace of Fi and 
F2 adiabatically when the Hamiltonian parameters change in time, allowing for an 
adiabatic transfer of the quantum state within the subspace. A qualitative measure of 
non-adiabatic effects is provided by the angle which is between S and its projection 
onto the subspace. 



Analysis of adiabatic transfer. 



9 



where Di = S ■ f j. When x ~ 0, the adiabatic following of 5" with the subspace is good. 
Following Oreg et al. [27], we can easily obtain x ^^i the dark state \ui{t)). The vectors 
5" and Fi for the above Hamiltonian are seen to be 

S= |^0,0,-sin(2^),0,0,0,cos2^,-^(l - Ssin^^)^ , 

f 1 = , ] = ( - gn, -nc, 0, 0, 0, 0, 

/gl + nl + 4AV3 ^ 

A, -A/VS 



The calculation is straightforward and we just quote the result for A = 0: Di = D2 = 
and X = 7r/2, thus making non-adiabatic effects rather strong. 

2.1. Fidelity of storage and retrieval 

We first consider the storage of a single photon from the cavity to the atomic memory. 
The initial state in this situation has ^lc{^) = ^0, and the atom-signal system is in the 
dark state (with zero eigenvalue), 

1^(0)) = cos^(0)|6,l) + sin^(0)|c,0) = |mi(0)), (27) 

where 

^(0) = tan-i (^) . (28) 

The signal pulse is to be stored by making the control field ^lc'{t) zero. A suitable form 
[20] for the control pulse is 

fi^(t) = fio [1 - tanh(rt)] . (29) 

The adiabatic evolution of the above state leads to just \ui(t)) at time t. After a lapse 
of time of the order 3/r, we expect the wavefunction to evolve to |c, 0) as fic(t) — )■ 
i.e., e{t) f . 

As a measure of any departure of our solution |^I^(t)) from the adiabatic answer, we 
compute the fidelity F(t) of the process given by 

F{t) = \{u,{t)\m)l (30) 

for different values of go and r. We also compute |y4o(t)P, |-Bi(t)p and |Co(t)P to portray 
the actual evolution of the state. 

For the retrieval of the photon from atomic memory, we need to increase the control 
field from zero to Qq. For this, we take the control pulse to be of the form [20] 

Qc{t) = ftots.nh{rt). (31) 

The initial state of the system is the dark state ( 127|) . now with 6(0) = |. Again, the 
fidelity fl30p of the process records the deviation from the adiabatic evolution. 
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3. Wavefunction formulation in the presence of dissipation 

There are standard ways of incorporating dissipation in quantum systems. A 
comprehensive account of those which are of particular use in quantum optics can be 
found in the textbook by Scully and Zubairy [28] . For a three-level system, the density 
matrix equations incorporating dissipation have been investigated in the literature in a 
somewhat different context [291 ED]. Here we adopt an approach which we believe to 
be quite transparent from a physical point of view, as supported by our results in the 
next section. We work directly with the wavefunction [31] , and this formulation can be 
regarded as an approximation to the full set of density matrix equations. We present a 
discussion of this approximation with regard to density-matrix treatment in Appendix 
A. Here we remark that our approach is similar in spirit to the formulation of Dalibard 
et al. and of Barchielli and Belavkin for a continuously measured system [321 l3T] . and 
has also been used in a similar context for A-systems in interaction with cavity fields 

In this approach [33], the unitary evolution of the system of a three- level atom 
interacting with the signal and control fields is interrupted by spontaneous decays. We 
assume that in the presence of spontaneous decay, which occurs over a negligible time, 
the system collapses to either level \h) or level |c) with equal probabilities. The decays 
occur in time according to a Poissonian distribution. To write down the wavefunction 
in this model, we use the following notation. U{t) denotes the evolution operator for 
the isolated system, which is computed in the previous section through the computation 
of Vra{t)s in fl20l) . The probability that a spontaneous decay occurs in the time interval 
dt is denoted by Tdt. The probability P{t) that a decay has not occurred for time t, 
after preparation of the system at t = 0, is e~^*. The operators that cause spontaneous 
decays to states \h) and |c) are, respectively, denoted by C,b and Cc- In writing down the 
wavefunction at time t, we have to include the possibilities of 0, 1, 2, . decays, 
with each of these weighted by the probability distribution mentioned above, and thus 

where Qis denote the possibility in which I spontaneous decays have occurred over the 
interval t. These are given as 

go = e^^*t/(t), 

= /* dti e-^(*-*^^f/(t - ti) r C e-^*^t/(ti), 
Jo 

ft fti ft2 fii — 1 

Qi= ■■■ dh dt2... dti e-r(*"*^) 

Jo Jo Jo Jo 

X u{t - ti) r c e-^(*i-*2)t/(ti - t2)T C . . . 

X e-^(*'-^"*')?7(t/_i - ti) V Ce-^'^U{ti), (33) 

where ( denotes either or (c- We shall also average over these stochastic histories 
by assuming the decays to be independent. In ( l32l) . Z{t) = w (\E'(t)|\E'(t)) is the 
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normalization of the wavef unction, which is necessitated as the evolution is no longer 
unitary. The summation over the series fl32l) in the present case is very easy, as the 
quantum evolution after the last collapse is what matters. The evolution from that 
state is either from level \b) or level |c). Suppose that the last Ith collapse occurred to 
1 6) level. Then 

Qi\^m= f dhe-^^'-''^VU{t-t,)\h) Pi{t,), (34) 

JO 

where Piiti) denotes the probability that / collapses have occurred in the interval to 
ti. This is given by 

Pi{h)=^^e-^'K (35) 



The summation over / now yields 

Y.Qi 1^(0)) = r / dt 
1=1 



xU{t-ti)\b). (36) 

Here we have used the result: 

oo 

Y.mi) = l-e~^'K (37) 
1=1 

One can write a similar expression if the last collapse occurred to level |c), by replacing 
\b) with |c). Combining these two possibilities with equal probabilities, we write the 
wavefunction as 

I^W) = -^^[e-''U{t)\m) + I [dh (l - e-^*^) 

X e"^(*-*i)f/(t-ti)(|6) + |c))]. (38) 

This is the final expression for the wavefunction in this model. 

To compute it, we again resolve it in terms of instantaneous eigenf unctions: 

\m{t)) =Y.Wu{t)e-'^>-'^'^''''\u,{t)). (39) 

k 

The coefficients Wk{t) can be straightforwardly expressed in terms of Vfc(t)s as 
W,{t) = ^ [e-^* V,{t) + \f^ dt, (e-^*^ - e-^*) 

i 

x((^(0)|6) + (M/(0)|c))], (40) 

where Xk{t) = jQXk{t')dt'. We solve these equations numerically and compute the 
fidelities (|5U]) and other relevant quantities in the storage and retrieval processes, as 
before, generalizing the coefficients V^s in flMl) - fl26l) to the above Wk{t)s. 
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t t t t 

Figure 2. Fidelity of the storage process versus time (in units of ^q^) with a signal 
Rabi frequency of go/^^o = 0.05: (a) T = 0, (b) T/fln = 0.1, (c) T/rJo = 0.5, and (d) 
r/ilo = 1. Each graph shows the effect of the variation of r: r/i7o = 0.1 (dashed), 
r/r^o = 0.2 (continuous), r/fio = 0.5 (dotted), and r/ilo = 0.8 (dot-dashed). 

4. Results and discussions 

We now present our results based on the equations developed in the last two sections. 
There are three relevant parameters which we take in the scaled forms of r/flo, Qq/VIq 
and r/^o- We have computed the variations of fidelities of the storage and retrieval 
processes with respect to all these parameters. The idea is to find the optimal parameters 
and develop an understanding of the dynamics of the transfer processes. A quantitative 
comparison with Berry's superadiabatic theory is presented in section 4.2 (see Table 1). 

4.1. Storage process 

For the storage of a single photon from the cavity to the atomic memory with the control 
field of the form fl29|l . figure [2] shows the fidelity f pOj) as a function of time (in units of 
fig ^) for a relative signal field Qq/VLq = 0.05 without dissipation (T = 0) at (a), and with 
dissipation at (b) T/Qq = 0.1, (c) T/Qq = 0.5, and (d) r/f2o = 1, for different rates 
of variation r (in units of Qq) = 0.1, 0.2, 0.5 and 0.8 of the control field in each case. 
Figures [3] and m show the same set of results for higher signal field strengths, go/^o = 0.1 
and 0.2, respectively. Note that the variation of the control field fl29|l ends as it drops 
to zero at t ~ 3/r, and indeed one finds that there is no change in fidelity values after 
this time. 

We first comment on the results without dissipation. In figure [2]^a), one sees that 
the fidelities saturate to higher values as r/Qo decreases, for example, F ~ 0.39 for 
t/Qq = 0.1 and F ~ 0.07 for r/Qo = 0.8. Thus fidelity of storage is better with a slow 
variation of the control field, as expected from the adiabatic theory. The same is true 
for other values of go/flo, as shown in figures [3]^a) andHJ^a). 

Next we examine how the fidelity depends on the signal strength. This variation 
is shown for a range of values of go/^o from 0.05 to 0.2, in figures El^ a) toH^a) for zero 
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Figure 3. The same as figure [2] but with a signal Rabi frequency of 50/^^0 =0.1. 
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Figure 4. The same as figure [2] but with a signal Rabi frequency of go/i^o — 0.2. 



dissipation. It is seen that the fidehty increases with an increase of Qo/^Iq, leading to 
F ~ 0.92 for r/fio = 0-1 ^'^ 9o/^o = 0.2. This can be understood easily as the non- 
adiabatic perturbation (proportional to 9) becomes large when flcit) 0, and it is 
larger, smaller the value of qq. 

We now examine the effect of dissipation. This is shown in figures |2]|1] in panels (b), 
(c) and (d) with increasing values of T/Qo. We see that with dissipation, the fidelities 
decrease, though the difference becomes marginal as F increases. The variation with 
respect to r/f2o follows the same trend of decreasing fidelities with increasing r/Qo, but 
is marginal for large F. Note that for a given system, F is a fixed parameter; however, 
here the dissipation rate is scaled by and hence its variation physically implies the 
inverse variation of the control laser power. 

To give a detailed picture of the evolution of our solution for the plots for 

the probability densities |ylo(t)P, \Bi(t)\'^ and |Co(t)P of finding the system in the states 
I a, 0), 1 6, 1) and |c, 0), respectively, are shown in figure [5] in case of 5^0/^0 = 0.1. On the 
left-hand side are the plots for no dissipation while on the right-hand side are those with 
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dissipation at T/Qq = 0.1. In storage, as Qq 0, the dark state approaches |c, 0). Thus 
|Co(t)P should be large at the end of the process. However, without dissipation, this is 
marginally fulfilled for the smallest rate r/Qo = 0.1, and in the presence of dissipation, 
it is worse. 

4-2. Retrieval process 

Next we turn to the retrieval of the photon from atomic memory with the control field 
of the form fl3T]) . Figure |6] shows the fidelity f pOj) of the retrieval process as a function of 
time (in units of ^q^) for a relative signal field go/^o = 0.05 without dissipation (F = 0) 
at (a), and with dissipation at (b) T/Qq = 0.1, (c) T/Qq = 0.5, and (d) T/Qq = 1, for 
different rates of variation r (in units of Qq) = 0.1, 0.2, 0.5 and 0.8 of the control field in 
each case. Figures [7] and [8] show the same set of results for higher signal field strengths, 
Qo/^o = 0-1 s-nd 0.2, respectively. We again note that the fidelity at a particular r/Qo 
saturates after a time t ~ 3/r when the control field variation saturates to Qq. 

As before, we first comment on the results without dissipation. We note from figure 
int^a) that the fidelities are much lower than that in the storage process in the range of 
values shown, even for rather small t/Qq. They do decrease with increasing t/Qq, but it 
is only a marginal effect. Thus the adiabaticity expectations are not quite fulfilled. This 
difference in the behavior of fidelity between storage and retrieval is easily understood 
by recognizing that for retrieval the non-adiabatic perturbation is the largest at the 
beginning of the procedure whereas for storage it is the largest toward the end of the 
procedure. 

As seen from figures [6](b)-(d), with dissipation the behavior is complex. For low 
dissipation rates (or high control power), the fidelity improves considerably but the 
behavior is non-monotonic. The best results are achieved with T/Qq = 0.1 and hereafter 
the results deteriorate with decreasing control power. Surprisingly, for retrieval we find 
that in the presence of dissipation, the general wisdom of adiabaticity is not followed; 
instead, the fidelity is better as r/Qo increases (except for very weak signal intensity, 
when Qo/^o ~ 0.01, not shown here). 

In figures [7| and [HI we mark the effect of the signal strength. It is seen from figures 
[MHlthat with an increase of the signal strength go/^o from 0.05 to 0.2, the fidelity of 
retrieval improves, as for storage. For each signal strength, there is an optimum value of 
F (or control power) at which the best fidelities are achieved. The behavior with respect 
to t/Qq shows the same unexpected trend. 

Again, to give a detailed picture of the evolution of our solution for {"^{t)), the plots 
for |Ao(t)P, |5i(t)|2, |Co(^)P, given in (El])- ([26]) with the coefficients V^s generalized to 
the Wk{t)s in f l^O]) . are shown for the retrieval process in figure [9]in case of Qo/^o = 0.1. 
On the left-hand side are the plots for no dissipation while on the right-hand side are 
those with dissipation at T/Qq = 0.1. For retrieval, as ~^ ^O) the dark state 
approaches |6, 1) (for VLq ^ go). So, |i?i(t)p should be large at the end of the process. 
However, as we find in figure [9] |i?i(t)p is quite small without dissipation, but improves 
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Figure 5. Plots of |Ao(i)P, |Co(t)P for the storage process, as a function of 

time (in units of fl^^) with a signal Rabi frequency of go/^o = 0.1 without dissipation, 
r = (left-hand side), and with dissipation at T/Qq = 0.1 (right-hand side). Each 
graph shows the effect of the variation of r: r/ilo = 0.1 (dashed), r/J7o = 0.2 
(continuous), r/ilo ~ 0.5 (dotted), and r/Og = 0.8 (dot-dashed). 
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t t t t 

Figure 6. Fidelity of the retrieval process versus time (in units of flg^) with a signal 
Rabi frequency of go/^o = 0.05: (a) T = 0, (b) T/no = 0.1, (c) T/fln = 0.5, and (d) 
r/f2o = 1- Each graph shows the effect of the variation of r: r/ilo = 0.1 (dashed), 
r/rto = 0.2 (continuous), r/fio = 0.5 (dotted), and r/flo = 0.8 (dot-dashed). 
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Figure 7. The same as figure [6] but with a signal Rabi frequency of go/^o =0.1. 
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Figure 8. The same as figure [H but with a signal Rabi frequency of go/^o = 0.2. 




Figure 9. Plots of |Ao(t)P, \Bi{t)\'^, | Co for the retrieval process, function of 
time (in units of fl^^) with a signal Rabi frequency of 170/^^0 = 0.1 without dissipation, 
r = (left-hand side), and with dissipation at T/Qq = 0.1 (right-hand side). Each 
graph shows the effect of the variation of r: r/ilo = 0.1 (dashed), r/fio = 0.2 
(continuous), r/ilo ~ 0.5 (dotted), and r/Og = 0.8 (dot-dashed). 



on inclusion of dissipation and it is better for a fast variation of the control field. 

It is finally noted that for slow rates of change of the control laser, the non-adiabatic 
corrections in the quantum evolution without dissipation can be estimated from the 
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Table 1. Comparison of our steady-state retrieval fidelities F for dissipationless cases 
(r ~ 0) with Berry's superadiabatic transfer results using (41). 
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superadiabatic transfer procedure given by Berry [2lj. This procedure has been apphed 
to a A-system with A = by Elk [3lj, using its equivalence to a two-level description 
[35] . Applying the above result for the probability Poo of transfer out of the dark state 
for the f2(7(t)-protocol f l3T]) for retrieval, we get 



The fidelity is then given by (1 — Poo)- We compare this result with our exact numerical 
computations in Table 1. 

As can be seen, our calculations are in good agreement with the superadiabatic 
theory, particularly for large Qo/^Iq. For small Qq and large r, the non-adiabatic 
perturbations are strong and we do not expect the superadiabatic formula to give reliable 
results. In general, for small r/Qo, the superadiabatic theory overestimates the fidelity, 
whereas at large t/Qq, it underestimates the fidelity. 



5. Conclusions 



In this paper, we have addressed the core problem of adiabatic transfer driven by 
a control laser for a single three-level atom confined to a high-finesse optical cavity. 
We have studied the problem numerically by using instantaneous eigenstates of the 
system, and obtained results with and without dissipation. The dissipation in our 
case is included by allowing for spontaneous decay of the uppermost level of the A- 
configuration. Without the inclusion of dissipation, our results concur with the adiabatic 
idea of increased fidelity for slower variation of the control field. The dependence on the 
rate r of variation is much weaker for the retrieval of the photon than for its storage, 
especially for weak cavity-atom coupling. The fidelities of both the processes are better 
for higher coupling at any rate of variation of the control field. 

In the presence of dissipation, the fidelity of storage still follows the adiabatic 
property; however, it gets worse with higher dissipation, particularly at slow variations 
of the control field. For a given cavity-atom coupling, the fidelity of the storage process 
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is better for higher ratios of control powers (Rabi frequency) to the spontaneous decay 
rate T. 

On the other hand, in the case of retrieval, the behavior with dissipation is non- 
monotonic. For an optimal control power Qq, whose value depends on the spontaneous 
decay rate F of the system concerned, the fidelity with dissipation is, in fact, better than 
that without dissipation, the best being at T/Qq ^ 0.1. Moreover, in the presence of 
dissipation, contrary to the adiabatic idea, retrieval fidelities are higher for faster rates 
of variation of the control field. 

In general, we conclude that an optimal dissipation aids the retrieval process but 
not the storage. Further, we find that the fidelity for storing is better, the stronger 
the control field and the slower the rate of its switching off. For the best fidelity of 
retrieval, however, there is an optimal power for the control laser dependent on the 
dissipation rate - the faster the switching rate, the better is the retrieval. Our results 
lend themselves to experimental tests. 

This quantum state transfer protocol, of course, has other practical limitations of 
implementation, e.g., residual atomic motion in the trap effectively reducing the atom- 
cavity coupling, imperfect preparation of the initial single-photon state in the cavity, 
etc. We have concentrated on the idealized model of a single atom trapped in a high- 
finesse cavity, which is relevant for studying the non-adiabatic issues discussed here in 
the presence of dissipation. 

For various quantum information processing protocols, one needs to either couple 
the three-level atom to external photon channels or deal with several atoms. The full 
quantum mechanics of adiabatic transfer in such situations is highly complicated due 
to the enlarged Hilbert space of quantum states as well as the more complex dissipative 
processes. We feel that our present analysis is a useful input for tackling such problems, 
and some work is in progress in this direction. 
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Appendix A. Density- matrix approach for dissipative systems 

We describe here briefly the dissipation model, used in section 3, for the evolution of 
the density matrix p{t). For a pure system, 

pit) = \J{t) p(0) Ut(t) 

= e-*«^* p(0), (A.l) 

where H^A = [H, A]. To describe the spontaneous decay, we introduce the notation 

pit) = c Pit) c^ (A.2) 

where can be Cb^ or Cc^- For example, 

Cb^p(t) = |fe)(&|. (A.3) 

Following the same approach as in the text, we include the possibilities of 
1,2,3,...,/,... decays, with each of these possibilities weighted by their probability 
distribution. Then we have 



pit) = Y. Piit)^ (A.4) 

1=0 



where 



and 



Po(t) = e-r* e-H'^* p(0), (A.5) 



Jo Jo Jo Jo 



lO JO JO Jo 
X rfti ^Xg-r(ti-t2) g-iHX(ti-fe) Y dt2C ■■■ 

X e-^*' e-'"'*' p(0). (A.6) 

As is well known [33], this result is equivalent to the following evolution equation 
dp 



dt ^ 



P 



= - i H V + F (Wp + pWt) + F WpW^ (A.7) 
where W = ( — 1. 

On the other hand, the wavefunction approach corresponds to the equation 

= H |vi;(t)) + zF(C - 1) |vi^(t)), (A.8) 
which in turn yields for the density matrix 

^ = -iH^ + F (Wp + pWt). (A.9) 
dt 

A comparison of flA.7p and flA.9p shows that the wavefunction approach agrees with the 
density matrix approach only up to flrst order in the 'no-decay' operator W. 
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